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Abstract 

The interpretation of experimental and numerical data describing 
off-equilibrium aging dynamics crucially depends on the connection be- 
tween spontaneous and induced fluctuations. The hypothesis that lin- 
ear response fluctuations are statistically subordinated to irreversible 
outbursts of energy, so-called quakes, leads to predictions for averages 
and fluctuations spectra of physical observables in reasonable agree- 
ment with experimental results [see e.g. Sibani et al., Phys. Rev. 
B74:224407, 2006]. Using simulational data from a simple but rep- 
resentative Ising model with plaquette interactions, direct statistical 
evidence supporting the hypothesis is presented and discussed in this 
work. A strict temporal correlation between quakes and intermittent 
magnetization fluctuations is demonstrated. The external magnetic 
fleld is shown to bias the pre-existent intermittent tails of the mag- 
netic fluctuation distribution, with little or no effect on the Gaussian 
part of the latter. Its impact on energy fluctuations is shown to be neg- 
ligible. Linear response is thus controlled by the quakes and inherits 
their temporal statistics. These findings provide a theoretical basis for 
analyzing intermittent linear response data from aging system in the 
same way as thermal energy fluctuations, which are far more difficult 
to measure. 
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1 Introduction 

How spontaneous fluctuations and linear response are related in off-equilibrium 
thermal dynamics is an open problem of considerable interest, as linear re- 
sponse measurements [H [21 [3l H] are the main in-road into the rich phe- 
nomenology of aging systems. Highlighting this issue are recent obser- 
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vations that aging dynamics is intermittent O O [71 |8]: rare but large 
fluctuations with an exponential size distribution punctuate much smaller 
equilibrium-like fluctuations with a Gaussian size distribution. Simula- 
tional [Hi nni [m [12] and experimental [H] evidence in difl'erent areas supports 
the idea that intermittent changes of magnetization and other observables 
are induced by, and hence statistically subordinated to, intermittent and 
irreversible outburst of heat, so called quakes. The same hypothesis leads 
to a widely discussed asymptotic logarithmic time re-parameterization of 
the aging dynamics (see e.g. Ref. [13j). To ascertain whether or not other 
mechanisms than subordination could produce these effects, direct statistical 
evidence is called for. 

In the following, statistical subordination and closely related issues are 
numerically investigated, using as a test-bed an Ising model, which is sim- 
ple and yet constitutes a bona fide instance of a complex aging system |11[ 
[T4l I15j. Intermittent changes in magnetization are shown to obey the same 
statistics as the quakes, and the magnetic field is shown to have negligible 
influence on the energy relaxation. Together, these two findings imply that 
aging dynamics is quake-driven. The fluctuation statistics is then investi- 
gated in detail, with the model properties in good agreement with previous 
investigations of intermittent heat flow [71 |Tl] and of magnetic linear re- 
sponse intermittency [8l[T2]. Quakes in this model are shown to be nearly 
uncorrelated events, which are described by a Poisson process whose average 
increases logarithmically in time. Additionally, a linear system-size depen- 
dence of the average number of quakes is found. This conflrms that quakes 
are spatially localized events, a property also found directly in this work via 
a real space analysis. The temperature dependence of the average number 
of quakes is very weak, except at the lowest temperatures. This hints to a 
hierarchical structure in the energy landscape of the thermalized domains 
spawning the quakes. 

2 Model and method 

In the model, Ising spins, cTj = ±1 are placed on a cubic lattice with 
periodic boundary conditions. They interact through the plaquette Hamil- 
tonian 



The first sum runs over the elementary plaquettes of the lattice, including 
for each the product of the four spins located at its corners. The second 



« = -E 




'Pijkl 



2 




t 8E 



Figure 1: (left): (colour on line) The average energy per spin, fig, is plotted 
versus the system age t, with the magnetic field switched on at times tw = 
200,400,600,800, 1000 and = 2000. The magnetic field has a quite smah 
effect on He- (right): The PDF of the energy changes 6E over a time 6t = 
100, with (stars) and without (circles) a magnetic field. Both PDFs feature 
a central Gaussian part of zero average and intermittent tails. The PDFs 
are nearly identical, except for the largest and rarest events. The simulation 
temperature is T = 1.5 for all plots. 

term describes the coupling of the total magnetization ^ cTj to an external 
magnetic field. As expressed by the Heaviside step function r]{tw — t), the 
field changes instantaneously at t = ty^ from zero to H > 0. Previous 
investigations of the model's properties in zero field show a low temperature 
aging regime |141 [TS] , during which energy leaves the system intermittently, 
and at a rate falling off as the inverse time |11| . 

The present simulations are all performed within the aging regime, i.e. in 
the temperature range 0.5 < T < 2.5, using the rejectionless Waiting Time 
Algorithm (WTM) [16]. The 'intrinsic' time unit of the WTM approximately 
corresponds to one Monte Carlo sweep. By choosing a high energy random 
configuration as initial state for low temperature isothermal simulations, an 
effectively instantaneous thermal quench is achieved. For each set of physical 
parameters, Probability Distribution Functions (PDFs) are collected over 
2000 independent runs, and other statistical data over 1000 independent 
runs. The symbol t stands for the time elapsed from the initial quench (and 
from the beginning of the simulations) . The symbol t^, is the time at which 
the field is switched on, while tobs = t — tw stands for the 'observation' 
time, during which data are collected. The value of the external field is 
set to H = 0.3 for t > tw The thermal energy is denoted by E, and 
the magnetization by M. The average energy per spin is denoted by ^e- 
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Figure 2: (colour online) (left): The PDF of the spontaneous magnetic 
fluctuations has a Gaussian central part and two symmetric intermittent 
wings. Data are sampled in the interval tyj^t^ + tobs- (right): The same 
quantity (outer graph, circles) when a field H = 0.3 turned on at tyj = 1000. 
The left intermittent wing is clearly reduced relative to the no-field case, and 
the right intermittent wing is correspondingly amplified. The inner, almost 
Gaussian shaped, PDF (stars) is the conditional PDF obtained by excluding 
the magnetic fluctuations which happen in unison with the quakes. 



The PDF of fluctuations in energy and magnetization are constructed using 
flnite time differences of E and M, taken over short time intervals of length 
5t < tobs- 

3 Results 

To qualify as a probe, an applied perturbation must not signiflcantly change 
the dynamics of the system investigated. Speciflcally for the present model, 
the average energy vs. time should be nearly unaffected. Secondly, fluctu- 
ation spectra should not undergo qualitative changes. Figure [H left panel, 
shows the average energy as a function of time, with six data sets correspond- 
ing to the perturbation switched on at tw values ranging from = 200 to 
tw = 2000. As expected, the field has little impact on the energy. The right 
panel shows, on a logarithmic scale, the PDF of the energy fluctuations in 
zero field (circles) and when a field is switched on at tyj = 1000 (stars). 
Gaussian energy fluctuations of zero average are flanked on the left by an 
intermittent tail, which carries the net heat flow out of the system. Again, 
only a minor difference is seen between the unperturbed and perturbed 
fluctuation spectra, and only for the largest and rarest of the fluctuations. 
Since the magnetic contribution to the average energy is negligible, quakes 
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Figure 3: (colour online) (a): The record signal, or Best So Far energy, 
for a simulation at T = 1.5 and no field. The full signal, which includes 
reversible fluctuations, is shown in the insert for a small sub-interval of the 
observation time, (h): active spins participating in any quakes occurring 
during the interval [10^, 1.1 x 10^] are shown in colour. All other spins are 
omitted. By definition, active spins change their absolute average value over 
five consecutive time units by exactly 2 (blue) or by 1.6 (green). Most spins 
are excluded as they either fluctuate with no average change or retain a fixed 
value during all 6Vs. 

the dissipation of the excess energy entrapped in the initial configuration, 
see Fig. [U in the same as in an unperturbed system would do [TT] . 

While a magnetic field induces a non-zero average magnetization, it does 
not change the overall structure of the spectra: The left panel of Figl2] shows 
the PDF of the spontaneous magnetic fluctuations occurring in the interval 
\tw,iw +iobs]- Intermittent wings symmetrically extend the central Gaussian 
part of the PDF. The outer curve (circles) in the right panel of the same 
figure depicts the PDF obtained when the field is turned on at = 1000. 
The positive intermittent tail is enhanced, the negative tail is reduced and 
the Gaussian part is not affected. Thus, the net average magnetization 
induced by the field, i.e. the linear response, arises through a biasing effect 
on the distribution of the spontaneous intermittent magnetic fluctuations. 

A key aging feature [3 [Tlj is that the excess energy trapped by the 
initial quench leaves the system through intermittent quakes. Intermittent 
magnetic fluctuations have a close temporal association to the quakes: The 
inner curve (stars) of Fig. [2] depicts a conditional PDF obtained by filtering 
out the magnetic fluctuations which occur simultaneously (i.e., in practice, 
either within the same or within the immediately following 5t ) with energy 
fluctuations of magnitude 5E < — 5. The filtering threshold utilized is near 
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Figure 4: (colour online) (a): The empirical cumulative distribution (stars) 
of the 'logarithmic waiting times' = In(tfc) — In(tfc-i) is plotted on a log 
scale. The full line is a least square fit to y oc exp(— ax), using the data 
points between the vertical lines, (h): The normalized correlation between 
Tfc and Tk+n (stars) is plotted versus n on a logarithmic scale. The line is a 
guide to the eye. We see that consecutive data points {n = 1) are only weakly 
correlated, and that the correlation decays exponentially for n > 1. This 
is in reasonable qualitative agreement with the theoretical approximation 
C{n) = 5nfi. 

the onset of the intermittent behavior of the heat flow PDF, as seen in 
Fig [TJ The filtering produces a nearly Gaussian PDF, demonstrating that 
quakes and intermittent magnetic fluctuations are synchronous events. In 
summary, the magnetic fluctuations which contribute to the linear response 
are subordinated to the quakes which dissipate the excess energy stored in 
the initial configuration. 

Visual inspection shows that energy traces of the p-spin model have a 
fluctuating part superimposed onto a monotonic step-wise decay, the latter 
given by the function rE{t) = miiiy^t E{t). This function is called record 
signal, or Best So Far (BSF) energy. The BSF energy is plotted in the 
main left panel of Fig. [3l The insert shows the full trace for a shorter 
interval of time. The right panel of Fig. [3] is a pictorial rendering of the real 
space positions of spins assessed to participate in one or more quakes during 
the time interval [10^,1.1 x 10^]. In the plot, the 'active' spins are shown 
in colour, and all others are omitted. By definition, spins are labeled as 
active if the time average over five consecutive time units of their orientation 
has changed by exactly 2 units (blue) or by 1.6 units (green). Note that 
the system-spanning spatial distribution of the involved spins integrates the 
effect of the many unrelated quakes occurring in the observation time. Any 
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Figure 5: (colour online) (a): The logarithmic rate of quakes a is plotted 
versus system size for T = 1.5. (b): The same quantity is plotted versus T, 
for N = 16^. The values of a are estimated from an exponential fit to the 
distribution of logarithmic waiting times, see Fig. [4^. 

single quake typically involves small clusters only. 

The quake statistics of this model should be amenable to a description 
as a Poisson process with average [IHl [HI [7] 

ni{t',t) = aln{t/t'). (2) 

Assuming that quakes can be treated as instantaneous events, which oc- 
cur at times ti, . . .tq . . ., a testable property equivalent to Eq. [21 is that 

the logarithmic differences (rather than linear differences) Tg =^ ln(tq/tq_i) 
are independent and identically distributed stochastic variables, sharing the 
exponential distribution 

Prob(rq > x) = exp(—ax). (3) 

A simple choice to evaluate the Tq is to identify the quake times tg with 
the steps of the BSF energy. Admittedly, this introduces an unwanted de- 
pendence on the data sampling frequency, when the latter is too high: the 
same transition from one metastable configuration to another may then be 
registered as multiple quakes. This leads to over-counting, and as discussed 
below, to spurious correlations appearing between successive quakes. 

From the empirical series of {Tg}q=i^2,... collected for each trajectory, the 
correlation function is estimated as C{k) = {TgTg^k)q — {Tq)q{Tq+k)q- The re- 
sult is then averaged over all trajectories. Systematic and statistical errors in 
the quake identification procedure lead to deviations from the theoretically 
expected Kronecker delta C{k)/C{0) = 6kfl- 
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In the left panel of Fig. H] the distribution of the Tq (dots) is well ap- 
proximated by an exponential decay over nearly three decades (line). For 
large values of the abscissa, deviations likely stem from a statistical under- 
sampling of rare events. For x ~ (the first three points are excluded 
from the fit) deviations indicate that successive intervals of equal or near 
equal duration appear more frequently than theoretically expected. This 
can arise when the interval between successive events is spuriously affected 
by the sampling time. The correlation function C(fc) of the logarithmic time 
differences, normalized to C(/c = 0) = 1, is shown in the right panel of the 
figure. After a drop to approximately 1/10 at A; = 1, the function tapers off 
exponentially with increasing k. The (modest) residual correlation confirms 
that successive events are occasionally mis-classified as quakes. 

A spatially extended glassy system with short range interactions is ex- 
pected to contain a number, say a, of independently relaxing thermalized 
domains, with a small and slowly growing characteristic size [20] '21] . Assum- 
ing that quakes originate independently within the domains, the observed 
fluctuations integrate the effect of a independent Poisson processes, whence 
a equals the logarithmic rate of quakes. Equivalently, its value can be es- 
timated from the exponential distribution of the logarithmic waiting times 
Tg, see Fig. IU as done in the present case. 

The number of domains, and hence a, must increase linearly with the 
system size, as demonstrated in the left panel of Figured) The temperature 
dependence (or lack thereof) of a is related to the geometrical properties of 
the configuration space, or energy landscape, of each domain. For record 
sized thermal fluctuations to elicit attractor changes, the energy landscape 
must be scale invariant [12]. As a consequence, changing the temperature 
should not change a at all. Note however that scale invariance cannot hold 
below a cut-off value where the granularity of the energy values makes itself 
felt. In the present model, the numerically smallest energy change following 
a single spin flip is 6E = ±4, i.e. unlike models with Gaussian quenched 
disorder, the granularity is important. The right panel of Figure [5] shows 
that a has a modest temperature dependence for 1 < T" < 2, i.e. for a major 
part of the range where aging behavior is observed. For lower temperatures, 
a clear T dependence is visible. 

4 Summary and conclusions 

Direct numerical evidence has been provided that intermittent magnetic fluc- 
tuations are statistically subordinated to a certain type of events, quakes. 
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which dissipate the excess energy trapped in the initial configuration. The 
external field does not alter the temporal statistics of the quakes and of the 
spontaneous magnetic fiuctuations. It only slightly biases the size distribu- 
tion of the latter. Therefore, the field can rightly be considered as a probe 
of the unperturbed off-equilibrium aging dynamics. In agreement with pre- 
vious investigations of other models [12] and experiments [8], the temporal 
statistics of intermittent energy and magnetization changes, is shown to be 
well described by Poisson process. 

Considering that aging dynamics is widely insensitive to details of the 
microscopic interactions, it seems reasonable to assume that the above find- 
ings are valid beyond the plaquette model itself. It should therefore be 
possible to analyze a wide range of intermittent linear response data from 
complex dynamical systems precisely as done for the intermittent heat flow 
data in Ref. 13 E]. 
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